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Abstract. 

In this paper a 3-phase Stefan problem solution method for ID semi-infinity alloy 
is developed. The problem is first solved for full enthalpy of the system and then 
the thermal difi'usivity has been eliminated from the divergence operator by Kirchoff 
transformation. Moreover, we introduce a similarity independent variable rj = x'^/t 
and original problem transforms to ordinary differential equations (ODE) for each 
phase separately. These ODEs have Dirichlet's type boundary conditions. 

Furthermore, we give an solution example of these equations for a simple case. 
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1. Problem and Full Enthalpy Equation. 

We formulate the problem for the temperature field T{x, r) (where x is a coordinate 
and T is the time) in following form. We have a ID semi-infinite liquid alloy array at the 
temperature Tmit > Ti, where Ti is the liquidus temperature of alloy. The alloy array is 
located at a; > 0. In the initial time r = the temperature at the point a; = instantly 
drops to Tout < Tg, where T, is the solidus temperature of alloy. Thereafter we have a 
consistent solidification process from left to right. We need to find T{x, r) and solidus 
Xs{t) and liquidus Xi{t) fronts locations. In any moment (except initial moment) there 
exist 3 phases: 

• Solid for < X < X,(r). 

• Solid-liquid (mushy zone) for Xs{r) < x < Xi{t) . 

• Liquid for x > Xiij). 

In this paper, we suppose p is constant and the same in each phase. Heat capacity C{T) 
and thermal conductivity k(T) are arbitrary single- valued functions of temperature T. 
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Mathematical the problem can be expressed as follows 

dT d\ 
pC{T)-^ = div {K{T)gradT) + P^-^-, (1) 

T{x, r = 0) = Tiniu (2) 

T(0, r) = Tout, r > 0, (3) 

T{oo, r) = Tinit, r > 0, (4) 

where L is latent heat, A is liquid fraction. We introduce the full enthalpy as 

T 

H{T) = pj C{t)dt + pL\{T) (5) 



then Eq. (P) can be rewritten in the form 
dH 

= div {K{T)gradT) . (6) 

With regard to thermodynamical identity 

gradH 

gradT = ^ , (7) 

dT 

we rewrite Eq. ® with initial and boundary conditions as 
OH 

= div {a{H)gradH) , (8) 
H{x, r = 0) = ffi^it = H{T,nit), (9) 
iJ(0, r) = Hout = H{Tout), r > (10) 
H{oo, r) = H.rnt = H{Tir,u), r > 0, (11) 

where a is thermal diffusivity. The temperature dependence of the thermal diffusivity 
is defined by 

«(^) = (12) 

dT 

The function (jH)) can be invert to T = T{H), therefore in (jH)) we write a(T) = 
a{T{H)) = a{H). It should be noted that the function a{H) has discontinuities at the 
Xg and Xi 

Oi{H)\H=H,-o = (^s, a{H)\^^jj^^f^ = ams, (13) 

a{H)\H=Hi-0 = Oirnh l/f=iy,+0 = (14) 

where we denote 

= H{Z), Hi = H{Ti). (15) 

Now, we use the Kirchoff transformation to eliminate a from below the divergence 
operator. We introduce the new function Z( X, t) as 

H 

Z{H) = [ a{h)dh. (16) 
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For Z{x, t) the original problem is given by 

Z{x, r = 0) = Zinit = Z{Hinit), (18) 
Z{0,t) = Zout = ZiHout), r>0 (19) 

Z{00, t) = Zinit = Z{Hinit), T>0, (20) 

where A is the Laplace operator, and a (if) = a{H{Z)) = P{Z), since we can invert 
Eq. (Uni) into H = H{Z). 

2. The Transformation of the Problem to ODE. 

Equation (|17|) can be transformed to a nonlinear ordinary differential equation (ODE). 
We introduce new independent similarity variable 

V = -, (21) 
r 

then we can write 



dZ _ 7] dZ dZ _ 1 dZ d'^Z _ 1 d'^Z 
dr 2t dr] ' dx -y/r dr] ' dx'^ r drf 

The Eq. fjl7|) can be written as 

d'^Z _ 7] dZ 
~d^ ~ ~2(3{Z)~d^' 
The initial and boundary conditions are 



(22) 



(23) 



Z{7] = 0) = Zout, Z{t] ^ Oo) = Zinif (24) 

Now, we suppose that the solidus/liquidus fronts positions are varying as 

Xsir) = hV^, Xiir) = kV^, (25) 

where kg and ki are constants. This dependence is the standard one for the ID Dirichlet 
heat transfer problem with phase change, which has many analytical and numerical 
solutions (see the references in [T]). Now we can divide the equation for each phase. Let 
(3{Z) be 

for [0,Z,), 

(3{Z) = {(3UZ) for Z e[Z,,Zi], (26) 

for Z G {Zi, oo). 
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where Z^^i = Z{Hsi). For each phase the equations have the uniform form 
d^Z V dZ 



d7]^ 2(3j{Z)dri 
The boundary conditions are 



i = s, m, I. 



(27) 



(28) 
(29) 
(30) 

Moreover, we have to define heat flux conditions on the interfaces between the phases. 
For Z these conditions are 



Z{r,= 


0) 


Zouti 


z{v = 




= Zs, solid phase {s) 


z{v = 






z{v = 


h) 


= Zi, mushy phase (m) 


ziv = 


h) 


= Zi, 


z{v^ 


oo) 


= Zinit, liquid phase (/) 



dZ 




dZ 


dx 




dx 


dZ 




dZ 


dx 


x=Xi-Q 


dx 



x=Xs+0 



x=Xi+Q 



pXoL— , 

dr 



(31) 
(32) 



where Aq - eutectic hquid fraction (for eutectic alloys). For the enthalpy we can write 



a. 



ami 



OH 
dx 
dH 

dx 



x=Xs-0 



x=Xi-0 



a. 



Oil 



dH 

IS T{ 

dx 
dH 

dx 



+ pXoL 



x=Xs+0 



dXs{T) 
dr 



x=Xi+0 



The unknown parameters kg and ki will be determined from either (jHT|) - 
3. Example. 



(33) 
(34) 

or 



The goal of this section is to give an example of the 3-phase Stefan problem solution. 
We will examine the problem with a constant thermal diffusivities in each phase 
tts 7^ ttm 7^ Analogous problems were studied in [TJ In this case we have 
(3j = aj for j = s, m, I. The general solution of the Eq. pTj) is given by 



Zj{ri) =Ci + C2y/7m]erf ^ 



(35) 



where Ci and C2 are arbitrary constants. Thus, it is to easy obtain expressions for the 
enthalpies {Hj = Zj/aj) in each phase 



H{x,t) = Hout + {Hs-H, 



out I 



f (2v/^) 
'f (2^) 



xe [0,x, 



(36) 



(Hi H,)erf {^J^) + H^erf {^X^) 








) 



(38) 



erfc ( — 5= ) 

H{x, t) = Hinit - (H.^it - Hi) 7 V^N ) X e (Xz, oo), 

. system for the determination of ks and ki can be written as 

( ) crf( ^ " ^P^ol^/c,, (39) 



The equation 
y/cTs^Hs - Hout)exp (-^) ^/c^{,Hi - H,)exp [■ 



Hi - Hs)exp (^-4^ j ^/^i{Hinit - Hi)exp (^-^) 
This system is easy to solve by numerical methods (see P and references therein). 



/a 



er 



(40) 
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